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ABSTRACT 

• In this paper we apply the jet-disk symbiosis model developed for Sgr A* to M81* - the 

Q> ' nucleus of the nearby galaxy M81. The model accurately predicts radio flux and size of M81* 

for the observed bolometric luminosity of the nuclear source with no major free parameter 
except for the inclination angle. We point out that the usually applied free, conical jet emission 
^ ^, model implies a longitudinal pressure gradient that must lead to a moderate acceleration of the 

■ jet along its flow direction. This, usually neglected, gradual acceleration naturally accounts for 

the inverted spectrum and the size/frequency relation of M81* and may be a general feature of 
radio cores. M81* is so far the best case for a radio-loud jet nature of the compact radio core in 
the nucleus of a nearby spiral galaxy. The fact that one can account for Sgr A* and M81* with 
, the same model by simply changing the accretion rate, strongly supports the jet-disk symbiosis 

' model as an explanation for the compact radio cores of galaxies in general. 
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tj^ , Quite a few nearby galaxies seem to have compact radio cores in their nuclei, prominent cases are 

' the Milky Way (Sgr A*), M31 and M81. Those radio cores resemble the cores of radio loud quasars, 

showing a very high brightness temperature and a flat to inverted radio spectrum that extends up to 
, submm wavelengths. Several models have been developed to explain those cores in the context of black hole 

' accretion: Melia (1992a&b) suggested a spherical accretion model for Sgr A* and, what he called, M31*. 

^ Falcke et al. (1993) and Falcke, Mannheim, & Biermann (1993, hereinafter FMB93) proposed an alternative 

model, where Sgr A* was explained as the core of a radio jet, fed by an underluminous, starving accretion 
disk (see also Falcke 1996a, for a review, and Falcke & Heinrich 1994 for M31*). This model evolved into 
the jet-disk symbiosis approach (Falcke & Biermann 1995, hereinafter FB95; Falcke 1996b) that unified the 
explanation for radio cores of quasars (Falcke, Malkan, & Biermann 1995, hereinafter FMB95), galactic jet 
sources (Falcke & Biermann 1996) and Sgr A* into a single picture. The basic idea was to postulate that 
black holes, jets and disks form closely coupled systems with little variations from one system to another 
except for the accretion rate. 

Recent VLBI (Bietenholz et al. 1996), multiwavelengths (Ho, Filipenko & Sargent 1996) and submm 
observations (Renter & Duschl 1996) have now shown that the radio core in M81 - which in analogy to Sgr 
A* and M31*, we will label M81* hereafter - is very similar to Sgr A*. This wealth of data now makes 
M81* an excellent laboratory to test the jet-disk symbiosis in detail. 
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2. Free Jet with Pressure Gradient 

One of the interesting findings of the VLBI observations (de Bruyn et al. 1976; Bartel et al. 1982; 
Bietenholz et al. 1996) is that the radio core of M81* is elongated and its size is intrinsically frequency 
dependent with r oc jy-o.8±o.05^ ^ frequency dependent size was one of the basic predictions of the jet 
model and excludes any homogenous, optically thin models (e.g. as in Duschl & Lesch 1994). However, 
the measured size-frequency relation is slightly shallower than in the Blandford & Konigl (1979, hereinafter 
BK79) free jet model adopted by FB95, and the radio spectrum of Sgr A* and M81* are both slightly 
inverted rather than flat as in the simple model. Here we show, that a self-consistent treatment of the BK79 
model implies a weak acceleration of the bulk jet flow due to its longitudinal pressure gradient which can 
naturally explain the observed size- and spectral indices in compact radio cores. 



2.1. Longitudinal Velocity Profile 

As in FB95, we shall describe the jet core within the framework of relativistic gasdynamics of a 
relativistic gas with adiabatic index T = 4/3. We consider only the supersonic regime and impose as 
boundary condition, that the jet expands freely with its initial soundspeed Ps without any lateral gradients 
behind the jet nozzle. This leads to the familiar conical jet with /An cx (BK79), where adiabatic 
losses due to lateral expansion need not be considered. For a confined jet scenario (e.g. Sanders 1983) at 
larger scales the adiabatic losses could, however, be quite severe with energy densities scaling as r~^^ . 

Nonetheless, even in a simple, free jet at least some work due to expansion will be done, because there 
is always a longitudinal pressure gradient, which will accelerate the jet along its axis. This acceleration is 
described by the z-component of the modified, relativistic Euler equation (e.g. Pomraning 1973, Eq. 9.171) 
in cylindrical coordinates where we set d/dr = and 8/ 86 = 

Here, uj — mp-nc^ + Uj + Pj ist the enthalpy density of the jet, [/j is the internal energy density, n is 
the particle density, and Pj = (F — l)C/j is the pressure in the jet (all in the local rest frame). For the 
'maximal jets' - the radiatively most cfhcient type of jet (FB95) - to be discussed here, we demand the 
equivalence of internal energy and kinetic (or better rest mass) energy Uj ~ nipnc^, hence uj = (1 + T)Uj 
and uj/n = (1 + r)mpC^ =const at the sonic point z — zq. In the free jet with conical shape the energy 
density evolves as Uj oc (7j/3j) and we do not consider any loss mechanisms other than adiabatic 

losses due to the longitudinal expansion. Using the relations mentioned above, the Euler equation becomes 




with ^ — (7j/3j/(r(r — l)/(r + l)))^^'". For ^ ~ 1 this is analogous to the Euler equation for the well known 
isothermal, pressure driven winds, with the wind speed replaced by the proper jet speed and the sound 
speed is constant and fixed by the required equivalence between internal and rest mass energy density at a 
value = V(r- l)/(r+ 1) - 0.4. 

The asymptotic solution of Eq. |^ for ^ =const and z — + cxd is 7j/3j oc 2\/ln z. If we ignore terms of 
the order ln7j/3j we could approximate the solution by 7j/3j ~ \/ rTT + 41n (z/zq)). For the following 
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calculations we will, however, use the exact, numerical solution to Eq. g ^^const), but the deviation from 
the approximate, asymptotic solution is rather small. 

2.2. Plasma Properties 

The basic ideas how to derive the synchrotron emissivity and basic properties of a jet in a coupled 
jet-disk system have been described extensively in FMB93 and FB95. The power Qj — qjAfdiskC^ — (Zj/i^disk 
of one jet cone is a fixed fraction of the disk luminosity, relativistic particles and magnetic field B are in 
equipartition within a factor fce+p - which we hereafter set to one, and the energy fluxes are conserved along 
the conical jet (for a moment we will forget the adiabatic losses). Mass conservation requires that the mass 
loss in the jet Mjet is smaller than the mass accretion rate in the disk Mdisk, thus q^^ = Afjot/Aifdisk < 1, 
while mass and energy of the jet are coupled by the relativistic Bernoulli equation 7j9m (l + /^s /(r ^ 1)) = 9j 
(FMB93). 

Here we will make use of the same logic and notation but with three changes with respect to FB95, 
namely (1) applying the velocity law Eq. H, (2) neglecting the energy contents in turbulence and (3) 
considering only a quasi monoenergetic energy distribution for the electrons at a Lorentz factor 70 (usually 
^ 100). The latter is indicated by the steep submm-IR cut-off of Sgr A* (Zylka et al. 1995) - and probably 
also M81* (Renter & Lesch 1996) ~ which is different from typical Blazar spectra and precludes an initial 
electron power law distribution (as in FMB93). 

The semi-opening angle of the jet \s — arcsin (7s/?s/7j/3j), and the magnetic field in the comoving 
frame of a maximal jet with the given sound speed, Ldisk — -^4i.5l0*^'^ erg/sec, gj/i = O.Sgo.Si a-nd 
^16 = z/lO^^cm becomes (see Eq. 19 in FB95) 

B = 0.6 G ^f}jL4i.5qo.5zT6 ■ (3) 

The number of relativistic electrons that are to be in equipartition with the magnetic field are a fraction 
Xc = n-c/n of the total particle number density, and the energy density ratio between relativistic protons 
and electrons is (/ip/c — 1)- From the energy equation 7jW7j/3jC7rr^ = ^j/i^disk we find that the characteristic 
Lorentz factor and electron density required to achieve equipartition are 

7e = m/ (ATmc^ip/cXc) = 344/ (/ip/cXo) (4) 
no = 45cm"^ /3jL4i.5go.5a;c^r6^ (5) 

Finally, to incorporate adiabatic losses we now have to make the following transitions with respect to the 
equations in FB94 

7j/3j ^ 7jA 7c ^ 7c • (7j/3j(^)/7j/?jM)'/' 

i3^i3.(7j/3j(z)/7j/3j(2o))'/^ (6) 

2.3. Synchrotron Emission 

The local spectrum of the jet will be Fi, oc i'^^^ between the synchrotron self-absorption frequency 
i'ssa(-z) and the characteristic frequency vdz) — 3ej^ suiacB{z)/A'!rmcC (here we will use an average pitch 
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angle Uc ~ 60°). For the maximal jet using Equations ||, | and |, we have in the observers frame 

v^{z) ^ 100 GHz V^J ^41.550.5/ (y\fPiTi^l^'p/cZiQj (7) 

with the Doppler factor T) = l/7j (1 — /3j cosi) for an inclination i of the jet axis to the line of sight. 

To find the local synchrotron self-absorption frequency v^sa we have to solve the equation 
T = ncisynciv) ■ "^r-J sini = 1 and transform into the observers frame. The jet radius is rj ~ (j)z and the 
synchrotron self-absorption cross section for a monoenergetic electron distribution is given by agync = 
4.9 • 10-"cm2 (S/G)2/37c"^/^(iy/GHz)-5/3^ yielding 



-'ssa 



.1/5.4/5 4/5 8/5 
3.3 GHz 2? ^^ ^^^41.5g0.5^o Mp/e ^ (8) 

7: (I?sini)3/5zi6 



If we now use the velocity profile Eq. |^ wc can invert Eq. ^ to find the size of the jet as a function of 
observed frequency. For the asymptotic regime z S> zqj J'c is approximated for a fixed inclination angle to 
within a few per-cent by a powerlaw in z, such that 



Zc = {zcfll s\ni){^p/^Xc) (90.5^41.5)^^^ ^^3.7^^10^3 (9) 

where the parameters are Z13.7 = zo/3 AU, z^io.s = 1^/22 GHz, ^ = (0.99,0.95,0.9,0.89,0.88) and 
Zc,o = (500, 1200, 1100, 900, 700) AU for i = (5°, 20°, 40°, 60°, 80°). 

The total spectrum of the jet is obtained by integrating the synchrotron cmissivity egync along the jet: 
(I'li) ^syncTT^dz, which Can again be approximated using powerlaws: 



^ 745mJy(qo.5i4i.5)^-^Vp%x^-^'^?3°?<f 

-337mJy(qo.5i4i.5)^-^Vp%^?-^^?3°?«, (10) 
= 247mJy(<zo.5i4i.5)^-"Vp%a;i-33z03i6i.o^i36 

-143mJy(go.5i4i.5)^-^Vp%:^^-''^?3'?^?o'l, (11) 

P — 1 nf; m T-irl'n T U. 40, ,.40 ^1.40 ^0.20, ,0.20 

ti, - iUbmJy((7o.5-t>4i.5j A^p/c^o ^13. 7^^10.3 

7lTnTir('n T U. 45, ,. 81 „1. 89-0.19, ,0.19 n 9\ 

- a mjy((70.5-t^41.5j A^p/o^^c ^13.7*^10.3; U^j 

for D = 3.25 Mpc, and i = (20°, 40°, 60°) respectively. This means that the resulting spectrum is no longer 
flat but tends to be inverted with a ~ 0.15 — 0.27 for i = 30° — 90°. The spectral index, like the size index, 
is a function of the inclination angle of the jet and tends to become flatter for smaller i. 



3. Application to M81* 

We can now apply this formalism to the situation in M81*, where we have five observed quantities, 
namely the time averaged spectral index a ~ 0.2 ± 0.2 (Renter & Lesch 1996), the flux Fiy(22GHz) ~ 120 
mJy, the size 2(22GHz) ~ 550 AU, the frequency dependence of the size z oc i/™ with m — — 0.8 ± 0.05 
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(Bietenholz et al. 1996), and the bolometric luminosity of the nuclear source (Ho, Filipenko, & Sargent 
1996) which is Ldisk lO"*^-^ erg/sec. 

Most of the model parameters are already well constrained. As the cosmic-ray ratio /ip/c and the 
fractional number of electrons Xe are almost interchangeable we keep the latter fixed at Xe = 1. The 
jet-disk ratio qyi was determined in FMB95 to be ^ 0.15 for a flat spectrum. To bring this model with an 
inverted spectrum and a peak at submm wavelengths to the same fluxscale one has to increase this value 
to qj/i = 0.5 where disk luminosity and jet have equal powers; here we do not intent to change this value. 
The size of the nozzle zq enters only weakly into the equations; if we assume that the high-frequency cut-off 
which is somewhere around 1000 GHz correspond to the size scale of the nozzle (in analogy to Sgr A* in 
Falcke 1996) we find by extrapolation that zq ^ 3 AU as an order of magnitude estimate. For a nozzle size 
of 10-100 gravitational radii this correponds to black hole masses of 0.3 — 3 • IO'^Mq which is consistent with 
dynamical estimates (Ho et al. 1996). 

This leaves us with two main input parameters: i and fJ.p/c. The cosmic-ray ratio determines the 
electron Lorentz factor and is well constrained; by definition fip/c > 1 and from the condition z^ssa < we 
here have ^ip/c ^ 5xc^^^ (Eqs. 0&|). Fi gure 1 shows the predicted sizes and fluxes of the model for various 
A*p/e and i, indicating that in M81* i ^ 30 — 40° and ^p/^ ~ 1.5 (FMB95 used ^p/^ = 2 for quasars). 
The predicted spectral index then is a ~ 0.15 — 0.19 and the predicted size index is —to ~ 0.9 — 0.92. 
Considering the simplicity of the model, those values are reasonably close to the observed ones. 

4. Discussion 

The jet-disk symbiosis emission model by FMB93 and FB95 which was initially developed to explain 
the radio core of Sgr A* and radio loud quasars can also explain the nuclear radio source in M81 in detail. 
Simply by scaling the accretion disk luminosity by several powers of ten to the observed value properly 
predicts flux, size, spectral index, and size index of the VLBI radio source. 

If one inserts M81* into the universal Ldisk/radio correlation presented in Falcke & Biermann (1996) it 
falls right onto the line that connects Sgr A* and radio loud quasars[|. This suggests that all these sources 
are powered by a very similar (jet/disk) engine. In this context it is quite interesting to note that M81* 
and Sgr A* are in spiral galaxies, yet appear radio loud in these diagrams. M81 also has recently shown 
double-peaked broad emission lines, a feature usually seen only in radio-loud galaxies (Bower et al. 1996). 

Like in quasars and Sgr A*, one has to choose the most efficient, equipartition jet model to explain 
M81*, with high internal energy and a jet power that equals the disk luminosity. This leaves one with 
no free parameters for M81* other than the jet inclination angle which according to the current model is 
around j ~ 30 — 40°. Some caution is, however, neccessary, as the PA of the VLBI component seems to 
change with frequency, something we have ignored here. It is not clear whether this indicates bending, 
helical motion, or an extrinsic effect and thus is difficult to interpret. Some of these effects could change the 
results slightly - but not to an order of magnitude. 

In this Letter we have considered for the first time the effects of the longitudinal pressure gradient in 
the BK79 and FB95 jet model. Without any sophisticated mechanisms, this gradient alone will already 
lead to a moderate acceleration of the jet to bulk Lorentz factors of 2-3. Consequently, the Lorentz factor 



^This was already shown in a preliminary way in Falcke (1994, Fig. 8.1). 
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will increase towards lower frequencies and for a fixed viewing angle outside the boosting cone the flux will 
become Doppler-dimmed with respect to higher frequencies. This can provide a natural explanation for the 
inverted spectrum seen in compact radio cores in general and for the size/frequency relation of M81*. Such 
a mildly accelerating jet may also be of interest for quasar and BL Lac radio cores or galactic jet sources. 
The biggest advantage, however, is that the jet velocity, which was previously a free parameter, now is fixed 
by a simple physical model. 

Another interesting point of the presented model is the usage of an initial, quasi monoenergetic 
electron distribution for which we find 7e ~ 200. This not only reduces another free parameter (the electron 
distribution index) but also can explain the submm-IR cut-off's seen in Sgr A* and M81*. If the jet in M81* 
does not have strong shoc;ks to re-accclcratc those electrons into the ususal power-law, the jet would have 
size-dependent high-frequency cut-offs and thus explain the apparent absence of extended VLBI components 
in M81*. 

This research was supported by NASA under grants NAGW-3268 and NAG8-1027. I thank an 
anonymous referee for helpful comments, and P. Biermann, W. Duschl, H. Lesch, P. Mezger and A. Wilson 
for numerous discussions. 
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